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Q I We use an exact holon and spinon Landau-liquid functional which describes the 

holon - spinon interactions at all scattering orders, to study correlation functions 

> 

OO ' of integrable multicomponent many-particle problems showing both linear and non- 



linear energy bands. Motivated by recent photoemission experiments, we consider 



Q^ ' specific cases when the dominant non-linear band terms are quadratic and apply our 

^' 

results to the evaluation of the half-filling ID Hubbard model one-particle spectral 

functions beyond conformal-field theory. 
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Recent experimental investigations on the one-particle spectral properties of quasi-one- 
dimensional (ID) insulators |1]] have confirmed that the low-energy physics of these materials 
is dominated by holons and spinous. The discovery of these elementary excitations followed 
from the diagonalization of the ID Hubbard model by Bethe-ansatz (BA). Theoretical 
studies of the one-particle spectral properties of this model are mainly numerical [|I],0,3. 
Analytical studies of the problem either focused on the limit of infinite on-site Coulomb 
interaction (when the charge - spin separation occurs at all energies) or to the metallic 
phase away from half filling, where conformal-field theory (CFT) has provided important 
information on correlation functions ^^- Unfortunately, CFT does not allow the evaluation 
of one-particle correlation functions for the half-filling Mott insulator and the study of these 
functions at finite values of U remains an open problem of great physical interest. 

In this Letter we combine a Landau-liquid functional which describes exactly the holon 
and spinon interactions at all scattering orders to study the effects of the non-linearity of the 
holon and spinon bands in the asymptotics of the ID Hubbard model one-particle correlation 
functions both at half filling and maximum magnetization. Our results apply directly to the 
study of the ID insulators and allow for the description of the unusual spectral properties 
detected in the low-dimensional materials |I|. 

We consider the ID Hubbard model in a magnetic field H and chemical potential n 
which describes N = Ni^ + N^ interacting electrons and can be written as H = T + U[D — 
^N] + Ea=c,sl^a'^S^- We have f = -tY.j,a[^]aCj+ia + h.c] the "kinetic energy" (t is the 
transfer integral), D = J2j %,!%,! measures the number of doubly occupied sites and /x^ = /i, 
/i^ = fioH (/io is the Bohr magneton) are the chemical potentials. Moreover, we define the 
charge and spin operators S'; = —^[Na — N], S^ = — |[iV| — iVj, N = J2a ^o-, N^ = Y.j %,ct 
where hj^a^ = Cj^Cja is the occupation number operator (c]^ and Cja are electron operators of 
spin projection a at site j = 1, ..., No). We define the charge density n = n^ + rii = N/Na 
and a spin density m = n^ — rii in the intervals < n < 1 and < m < n, respectively, 
where n^ = N^j/Na- The relevant Fermi momenta are kpa = Trno- and kp = TTn/2. We use 
units such that a = h = 1, where a is the lattice constant. 



In the operational BA representation of Ref. [0 the holons and spinons are the holes 
of the c and s anticommuting pseudoparticles P, respectively. In this basis the Slater- 
determinant levels of a pseudoparticles, with a = c,s, are exact energy eigenstates of the 
many-electron problem and the total momentum is additive in the pseudoparticle momenta 
q. The Hamiltonian contains pseudoparticle interaction terms, yet integrability implies that 
all the interactions are of zero-momentum forward-scattering type. The exact many-electron 
ground state (GS) is in this representation described by the a-pseudoparticle distribution, 
K{q) = 0(#J^ - q) and iVO(g) = e(g - gS^J^) with < g < qi+^^ and qi'^^ < g < 0, 
respectively, where qp^ = ±qFa + 0{l/Na) and qpa = T^Na/Na- Here Nc = N and A^^ = A^j, 
and hence qpc = '^kp and qps = kp^. The pseudo-Brillouin zones are q^^^ = ±qa + 0(1 /No) 
and qa = nN^/Na where A"* = Na and A"* = N^ and hence gc = tt and g^ = kp^. Note that 
at half filling, A'^'* = N*- Nc = Q (zero magnetization, A'^'* = A"; - A'^ = 0) and the GS has 
no holons (spinons). Importantly, Aq, and A"* are always conserving numbers. The energy 
eigenstates only involve the momentum distribution Na{q) and for states differing from the 
GS by a small density of excited a pseudoparticles the excitation-energy spectrum can be 
expanded in the momentum deviations, AA'Q,(g) = Na{q) — N^{q), what leads to a Landau 
functional |^, AE^ = J^'iZi AEi, which up to third order in i reads 

AEp = Y.AN^{q)e^{q) + ^ ^ ^ AAr,(g)AAr,,(g')^ ^,„,(g, g') + 

a,q ^ a,a' q,q' 

+ A E Y. ^NM^N^'{q)^N^"{q")]:9a,c.',AqA'.q") + h.o.. (1) 

a, a', a" q,q' ,q" 

All the coefficients of the functional (|I|) can be exactly derived from the BA equations [§ , 
yet for orders i > 2 the calculations are very lengthy. So far only the i = 1 band ea{q) and i = 
2 fa,a' function were derived |^. /^From the band expressions the velocity fa(g) = dea{q)/dq 
and the function aa{q) = dva{q)/dq follow. The fa^a' function expression only involves 
the velocity and the two-pseudoparticle phase shift ^a,a' {q, i') • This phase shift plays 
a central role in our critical theory - it is an interaction-dependent parameter associated 
with the zero-momentum forward-scattering collision of the a and a' pseudoparticles of 
momentum g and g', respectively. For instance, the phase-shift combinations ^^^.z = 5a,a' + 



J2j=±iUy^aa'i(lFa, JQPa') with i = 0,1 fully determined the a anomalous dimensions of the 
generahzed critical theory, as we discuss below. 

When both the densities of a pseudoparticles, Ua, and of a pseudo holes, n'^, are finite 
the dominant €a{q) band terms are linear in the momentum for low-energy excitations. We 
emphasize that in this case only the terms of order i = 1 and i = 2 oi the functional (0) are 
relevant for the low-energy physics. Upon linearization of the ea{q) bands and replacement of 
fa,a'{(l, q') by fa,a'{(lFa, ^QPa'), the spcctrum (|l]) leads to the CFT spectrum of Ref. p. This 
reveals that our basis is the natural representation for the critical-point operator algebra. 
For these densities CFT provides asymptotic correlation functions expressions [^. 

The functional (|l|) is universal for a wide class of integrable multicomponent models. As 
we illustrate below for the ID Hubbard model, occurs often in these systems that as the 
density ria or n'^ becomes small the dominant low-energy a band term is non linear and of 
the general form ea{q) ~ ilkl —Qa\-'/j^-M.a, where Qa = when na — *> and Qa = Qa when 
ra^ — > 0, M-a is a positive generalized mass, and j is a positive integer, j = 1, 2, 3, 4, ...., cxd. 
Thus, j characterizes the curvature of the band at the pseudo-Fermi points, that is, j = 1 
is linear (the usual case described by CFT), j = 2 is quadratic and so on (j = oo describes 
the dispersionless band). We use the notation {a,j} for a a band of type j. The c and s 
bands of the present model are represented in Figs. 7 and 8, respectively, of Ref. 0. For 
U > the c band is quadratic around q = ±7r and q = and at m = and U > the 
s band is only quadratic around q = 0. In the limit oi U ^ oo its width vanishes and 
it becomes flat, whereas the width of the c band is always At. The c and s pseudo- Fermi 
surfaces are at g = ±2kF and q = ^kpi, respectively, and at low energy the bands can be 
linearized around them outside the quadratic or flat regions. Therefore, there occur band 
changes {c, 1} -^ {c, 2} as ric -^ and n^ ^ when 2A;^ -^ and 2kp -^ n, respectively, 
and {s, 1} —>■ {s, 2} a.s rig —>■ when kp^ —>■ 0, and {s, 1} —>■ {s, oo} as f/ ^ oo at m = and 
the s band becomes flat (es(g) = implies {s, oo}). For each {a,j} branch the low-energy 
critical spectrum has contributions from terms up to order i = j + 1 oi the functional (|T]). 

In this Letter we study the effects of non-linearity on the asymptotics of correlation 



functions. In particular, following recent photoemission data |]l|, we consider the one a 
electron (or hole) spectral function A^ik.uj), momentum distribution N„{k) = I]^ v4o-(A;, cu), 
and density of states Dcr{uj) = J2k^cr{k,uj). There are no previous analytical results for 
these functions at finite values of U and (I) for the n = 1 Mott-Hubbard insulator and (II) 
the fully polarized ferromagnetic {m = n) initial GS's. At n = 1 and m = there are 
neither holons nor spinous in the initial A^-GS and the final A^ — 1 states have one holon 
and one spinon. Importantly, the numerical energy dispersions of Ref . [|l| , which agree with 
the phtoemission data, can be described by our holon and spinon bands ea(g). At large U 
the right (or left) holon line, which we denote by uj^{k) [or uj^{k)]^ of Fig. 3 of the first 
paper of Ref. |jl[ is in the thermodynamic limit slightly uo shifted and refers to c<j > (or 
cj < 0) and reads uj^ik) = 2t + ec(vr/2 + k) [or oj^ik) = 2t + ec(7r/2 — k)]. It corresponds 
in our representation to the creation of the spinon at g = — 7r/2 (or q = 7r/2) and the holon 
in the domain |g| G [tt/2,tt] (or |g| G [0,7r/2]). We emphasize that the final (A^ — 1)- 
GS corresponds to k = 7r/2 = kp and (due to the slight uj shift) to uj = 2t. The right 
holon line is quadratic in u around the GS energy because of the j = 2 character of the c 
band at g ~ ±7r. To illustrate the physical importance of our non-linear critical theory, we 
evaluate below the u dependence of the uj^{k) peak which in the A"a — > oo limit exists at 
k = kp = vr/2, between the k = Svr/T and k = An/T peaks of the figure. [Below we measure 
u from the GS, i.e. uj^{k) = ec(7r/2 + k).] We also find that the n = 1 non-linear effects 
lead to a non-Luttinger-liquid divergent behavior for D^^u). 

The low-energy critical theory with both {a, 1} and j > 1 {<y,j} bands (for our model 
the total number of bands is two) requires the use of terms up to the order i = j + 1 
in the functional (|l]). For i = 3 the calculations are lengthy and are technically similar 
to the the ones presented in the Appendix of Ref. [^ for order i = 2. Here they involve 
substitution of the third-order deviation expansions in the BA equations of Ref. [§[. This 
leads to ga,a',a"{q, q', q") = ga,a',a"{q, q', q") + ga,a',a"{q', q", q) + ga",a,a'{q", q, q') where 

ga,a',a"{q, q' , q") = 4:{7lf{aa{q)^a,a'{q, q')^a,a"{q, q") 



j=±l a'" "^ 

+ Va{q)[^a,a"{q,q")Y^»,a'{q,q') + Aa,a',a"{q,q\q")] 

+ XI ^Va"'^a"'AJqF(x"',q)Aa"',(x',a"UqFa"',q\q")}, (2) 

i=±l a'" 



Aa,a',a"{q,q',q") = ^a,a' {q, q')-r^o,,a" {q, q") 

dq 



d ,.$«/,«(g',g)<l'«/,„"(g',g") , $„",„(?", g)$a",a'(9", ?' 



+ 2vrp„ {r^\q)) f {[ ""-" ^'^ ' ^g;"");/ '^ ^ + 



^g 2vrp„,(r2)(g')) 27rp„» (r2)(g")) 

jt±i^ 2vrpa/// (ro,///) 

Wq = Vaiqpa), Oq = ciaiqFa)-, the functioii 2Tcpa{r) is such that l/27rpa (r(j°)(g)j = 
drj^\q)/dq, rf\q) = K^^\q), rf\q) = S^^\q), K^^\q) and S^^\q) are the rapidity functions 
obtained by solving Eqs. (Al) and (A2) of Ref. g| for N'^{q), and Tq, = 'r^°''((3'Fa)- 

The Landau character of the spectrum (^ imphes that the only small parameter which 
defines the range of validity of the corresponding low-energy critical theory is the density 
of excited pseudoparticles. The associate low-energy processes can produce changes in the 
numbers A^'^^ of right (i = +1) and left {t = —1) pseudoparticles (with l = sgn{q)l = 
±1), which are conserved independently. Equivalently, they can produce changes in the 
charge and current numbers Na = J^l^m and Jq = (1/2) J2l '•^qm respectively. The above 
processes can be classified into three types, (i) GS - GS transitions associated with variations 
AA^Q, = Na — N^ which only change qpa, (h) finite momentum K = J^a'^a'^qFa processes 
associated with variations V^ = Ja — Jay and (iii) a number N^^ = 0, 1, 2, ... of elementary 
particle - hole processes around the Fermi point qpa- In the absence of transitions (i) we 
define qpa relatively to the initial GS and otherwise relatively to the final GS. Independently 
of the form of the band ea{q), at low energy the above processes lead to a momentum of 
the form K = ko + J2a ^Pa where AP^ = ^[^NaVa + Y,l '■^af] and in the present case, 
^0 = HaT^a^qFa- Let US deuotc by X^s, Ha^ and Y,a the summations over linear branches, 
quadratic (or other j > 2 non-linear) branches, and all types of branches, respectively. 
The {a, 2} band reads ea{q) = ^\\q\ — QQp/2"^a5 with Aia = ^a ^^"^ "^a — V'^a(O) ^^^ 



m* = l/aa{qa) for pseudoparticles and pseudoholes, respectively. These masses play the 
same role as the effective mass of Ref. [|1^]. For the Hubbard model m* and m* are given 
by Eqs. (47) and (48), respectively, of Ref. p. Introducing in Eqs. (|l])-(il) the {q^)J} 
band expressions (with j = 1 or j = 2) and the low-energy deviations associated with all 
low-energy processes and performing the q summations leads to the critical spectrum 

^E = ^E^.Af + i:UK?] + E{(^«)'^«)' + (^)'% + ^^«Arf } , (4) 

L/ ^ „ „ 07T L/ Tn„ Li 

where Af = T.^=o^{K? + iVf , K = E.'C'^, A^ = E.'^L'^a', v^ = v^ea with 
Va = qpa/ma and Va = qpa/m*^ (or Va = [qa - qpal/rfia and f„ = [q^ - qFa]/m*J, ^° = 
^ = ^ = -^j^, and £a is independent of -Uq and contains no pseudoparticle interaction 
parameters. Note that qpa = T^^N^/Na (or q^ — qpa = —T^^Na/Na) in the limit when 
n^ = (or n^=0) for the initial GS. 

When all a bands are linear, the critical-theory expressions involve the anomalous di- 
mensions, 2A^ = Y.i=Q,i{.^if + i^NJ)o. + 2NP^^^ (equal to (Ei=o,i '^"^if + '^Np^^, when there 
is dual symmetry) and the asymptotics of correlation functions is given by CFT [^ and reads 
Xtf(x, t) oc na,Je~*'^"^/(a; — '•fat)^^"]. On the other hand, the critical theory associated with 
the low-energy spectrum (^ is characterized by two remarkable properties: (I) in spite of 
its non-linear velocity terms, the same above anomalous dimensions are obtained from the 
following equation, 2A^ = ^[Sa^ad^E /dva + 5a,adA.E/dva + lAPa]. Since their expres- 
sions are the same as in the linear regime, their values are given by taking the corresponding 
limits of the densities in the general expressions - in the present model the values of the 
dimensions 2A^ (or 2A^) are in the quadratic case given by the limit n -^ 1 (or m ^ n) oi 
the corresponding dimensions of the linear regime; and (II) a.s n^ —>■ or n'^ —>■ there is an 
adiahatic symmetry |11| which implies that ^°„/ , ^\i^ -^ 5a,a' in these limits. (However, if 
the a' band is linear, E,^,^ and E,\^i remain in general interaction dependent.) This symmetry 
leads to ^° -^ 1, i.e. the pseudoparticle interactions which renormalize the quadratic mass 
m* and associate velocity Va vanish (and rha = fn'^). For instance, in the present model 
Ca ' Cc -^ ^c,a asn^l {i.e., n'^ -> 0) and 4°^ , ^l^ -^ 6s,a as m ^ n {i.e., n^ -> 0). 



These symmetries imply that the terms of the critical spectrum (Q) which contain no lin- 
ear velocity terms are of non-interacting pseudoparticle character. Following this remarkable 
property we find the following expressions for the asymptotics of correlation functions 

oc n l/[(-^Wt)'^^' (v'2t/m*)S<2A^], x = 0. (5) 

These x and t dependences can be understood in the following way. In the limit of low energy 
each a excitation branch corresponds to an independent momentum-energy tensor compo- 
nent and to one independent Minkowski space with light velocity Va- For bands with both 
finite Ua and n^ densities the velocity v^ is also finite and the metric is Lorentzian. However, 
for vanishing small values of n^ or n^ the metric becomes Galilean. In the Minkowshian 
case the asymptotic of correlation functions involves the variables {x ± Vat) associated with 
Lorentz transformations. On the other hand, in the case of Galilean symmetry space x and 
time t are transformed independently. This is consistent with the two asymptotic-expression 
regimes involving either x or t. Moreover, the only combination of the time t and mass ?7i* 



with dimensions of x is const x Jt/m*^. Note that the changes in the asymptotics only 
concern the metric whereas the a anomalous dimensions, whose values depend on [//t, n, 
and 771, remain the same. Both the CFT asymptotic correlation function expressions and 
expressions (|^) are limiting cases valid for finite and vanishing, respectively, values of n^ or 
n^. As Ua or n^ is gently increased, we come into a small-density {a, 2} -^ {a, 1} band 
transition regime which is not described by these asymptotic expressions. 

Our theory does not describe the case when all bands are of j > 1 type. When some of 
the bands are of j > 2 type we find expression (|^) for x ^ {t/AiaY^-', whereas the x = 
expression becomes Xisi^^'t) oc Ila.a'/' l/ll—'-'^a'^)^ "' (t/Aia) ^ ]• Here j is meant to be 
a function of a, i.e. different a bands may have different j > 1 values. 

Fourier transforms of the above asymptotic expansions provide correlation-function ex- 
pressions for values of momentum k close to ko and low values of energy u measured from the 
initial GS energy. For 171 = and both n = 1 and small finite densities of holes 6 = {1 — n) 

8 



and low negative (positive) values of uj and/or values of k close to kp our generalized theory 
leads to A^jikp^uj) oc |u;|~'^/^ and N^(k) oc \k — kp\^^^ for the particle (hole) spectral function 
and momentum distribution. (At ra = 1 we define the ground-state energy at the bottom 
and the top of the Mott-Hubbard gap [Q for particles and holes, respectively.) We emphasize 
that Af,.{kF,uj) oc |c<j|^'^/^ is (for particles) the above n = 1 peak of the figure. 

On the other hand, the particle (hole) density of states is given by D„{uj) oc |u;|^^/"'^^ 
and D^{uj) oc \uj\^'^ for n = 1 and small finite 5 values, respectively. However, the latter 
expression (also predicted by CFT) is restricted to frequencies \uj\ < Ec = ^(x)^ 32m* - 
In the limit of n — > 1 this domain shrinks to a single point and the spectral function 
diverges as D{uj) oc IcjI"^/^^. In the m = n case we consider creation of one | electron 
and find for both m = n and for small finite densities n^, A^^kpi^uj) oc |ct;|^^+2[i^''ol and 
N^{k) oc |A; — A;^||2[^^''o] , where t]q = (2/7r) tan^^([4tsin(7r?2)]/?7) and n < 1, whereas the 
density of states is given by Di{uj) oc jco'l^^+^Ii"'"'] a^d D^^u) oc Iculal^-w] for m = n 
and small finite values of n^, respectively. The latter expression is restricted to energies 
\uj\ < Eg = <^(^)^^- Again, in the limit of n^ -^ this domain shrinks to a point 
and the spectral function behaves as -Dj.(co') oc |co'|^2 + 2[i^''o] . These Di(uj) expressions are 
not valid for U ^ oo because then the bands are of {c, 1} and {s, oo} type and instead 
Di{uj) oc |c<j|~"2. Finally, concerning the comparison of our results with previously obtained 
m = and U —>■ oo expressions, while our theory does not apply to the n = 1 case of bands 
{c, 2} and {s, oo}, it provides expressions for the {c, 1} and {s, oo} case which corresponds 
to finite values of n and S. Importantly, our general expressions lead to the same results as 
Ref. 0, with the density of states given by Dtj{uj) oc Itul^^/^ and D^{uj) oc \uj\^/^ for U ^ oo 
and small finite values of 4:t/U, respectively [|I2|. Our non-linear critical theory results are 
expected to shed new light on the unusual properties of quasi-lD materials. 
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